jsMath

Uvjetovanost Hilbertove matrice

def hilbert_matrix(n): return matrix([[1/i for i in [j..j+(n-1)]] for j in [1..n]]) 
       
def hilbert_cond(n): A = hilbert_matrix(n) return A.norm()*A.inverse().norm() 
       
hilbert_cond(5) 
       
476607.250243
hilbert_matrix(30).inverse().norm() 
       
2.12827697168e+43